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The aim of this paper is to characterise the star formation 
activity in the poorly studied embedded cluster Serpens/G3- 
G6, located ~ 45' (3 pc) to the south of the Serpens Cloud Core, 
and to determine the luminosity and mass functions of its pop- 
ulation of Young Stellar Objects (YSOs). Multi-wavelength 
broadband photometry was obtained to sample the near and 
mid-IR spectral energy distributions to separate YSOs from 
field stars and classify the YSO evolutionary stage. ISOCAM 
mapping in the two filters LW2 (5-8.5 /mi) and LW3 (12- 
18 yum) of a 19' x 16' field was combined with JHKj data 
from 2MASS, K s data from Arnica/NOT, and L' data from 
SIRCA/NOT. Continuum emission at 1.3 mm (IRAM) and 
3.6 cm (VLA) was mapped to study the cloud structure and 
the coldest/youngest sources. Deep narrow band imaging at the 
2.12 fmx S(l) line of H 2 from NOTCam/NOT was obtained to 
search for signs of bipolar outflows. We have strong evidence 
for a stellar population of 3 1 Class II sources, 5 flat-spectrum 
sources, 5 Class I sources, and two Class sources. Our method 
does not sample the Class III sources. The cloud is composed 
of two main dense clumps aligned along a ridge over ~ 0.5 pc 
,— ' | plus a starless core coinciding with absorption features seen in 

■ the ISOCAM maps. We find two S-shaped bipolar collimated 
I flows embedded in the NE clump, and propose the two driv- 

ing sources to be a Class candidate (MMS3) and a double 
ON , Class I (MMS2). For the Class II population we find a best age 
of ~ 2 Myr and compatibility with recent Initial Mass Functions 
' (IMFs) by comparing the observed Class II luminosity function 

■ (LF), which is complete to 0.08 L Q , to various model LFs with 
\Q " different star formation scenarios and input IMFs. 

o \ 
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1. Introduction 

The Serpens Cauda Clouds are among the dark est regions in a 
comp lex of clouds called the Aquila Rift JPame & Thaddeusl 
Il985h . The young embedded cluster in the Serpens Cloud 
Core has been well studied at most wavelengths, see 
lEiroa. Diupvik & Casalil Efiwk for a recent review. About 45' 
to the south th e re is a lso a region with active star formation. 
ICohen & Kuhl dl979h found 4 op tical T Tauri stars in a 30" 
field and named them Ser/G3- G6. IClarkl dl990lll99lh mapped 
the area surrounding Ser/G3- G6 in the ammonia 1,1 emission 
line and found two NH3 cores, one on each side of the opti- 
cally visible stellar group : Ser/G3- 6NE and Ser/G3- G6SW. 
IZiener & EisloffeJ i 19991) surveyed 3.14 square degrees in the 
optical lines of [S II] and found the Herbig-Haro object HH 476 
close to the Ser/G3-G6SW core. T he energy source of HH 476 
was identifie d bvlWu et alJ J2002h to be I RAS18265+0028. At 
this position Persi, Palagi, & Felli ( 1994) found a H2O maser. 
Recently, Ha rvev et alJ ( 120061) found 24 YSO candidates in a 
12' x 12' area in this region (which they refer to as Cluster B) 
from a mid-IR survey with Spitzer. 

When preparing the ISOCAM star formation survey . 1 this 
region was defined based on the ammonia maps of Clark 
( 1990) and named Serpens NH3. The field is relatively opaque 
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at optical wavelen gths and has a peak brightness at 60 pm 
JZhanget alJl988h . 

We will use the distance o f 225 + 55 pc to the Ser pens 
Cauda Clouds according to IStraizvs. Cernis. & Bartasiute I 
d2003l) . keeping in mind the caveat of the same authors that 
the cloud complex is possibly 80 pc deep. 

Preli minary ISOCAM r esult s from this reg ion were pre- 
sented in iKaas et alJ dl999l) and iKaasI Jl999b). In this paper 
we present the full results of a 17' x 19' ISOCAM survey in 
LW2 (6.7 pm) and LW3 (14.3 pm), cross-correlated with the 
2MASS point source catalogue. In addition we present deep 
follow-up H2 line (2. 122 pm), K$ band (2.15 pm), and L' band 
(3.8 pm) imaging from the Nordic Optical Telescope, as well 
as IRAM 1.3 mm continuum and VLA 3.6 cm mapping of the 
most active part of the region. 

We describe the various observation sets with reductions 
and results in Sect. |2] the use of IR excesses to extract YSO 
candidates in Sect. [3] the classification of the YSO population 
in Sect. [4] the luminosity and mass function in Sect. [5] and the 
cloud structure and the outflow sources in Sect. [6] 

2. Observations and reductions 

2.1. ISOCAM 

This paper is based on observations obtained with ISOCAM, 
the camera aboard t he Infrared Space Ob servatory (ISO; see 
lKessleretaTlll996l ICesarskv et~al1ll996l) . and is part of the 
LNORDH.SURVEY_l star formation survey conducted in the 
two broa dband filters LW2 (5-8.5 u rn) and LW3(12-18 ^m), 
see, e.g., Kaas & Bontemps (2000) for a review. In the re- 
gion called Serp-NH3 about 0.09 square degrees were sur- 
veyed in LW2 and LW3 in 1996. The pixel field of view 
(PFOV) was 6" and the unit integration time 0.28 seconds. 
The region was mapped by rastering along right ascension with 
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about half a frame (90") overlap in a and 24" overlap in 5, 
producing a 12 by 6 mosaic. Each position in the sky was 
observed for about 15 seconds. Image reductions were per- 
formed using CIA 2 and a set of our own p rograms for point 
source detection and photometry. We refer to lKaas et al. ( 2004) 
for a detailed description of the data reduction. The fluxes 
in ADU/s in LW2 and LW3 are converted to mJy through 
the relations 2.32 and 1.96 ADU/gain/s/ mJy for LW2 and 
LW3, respectively JBlommaert et alJl2000l) . These are strictly 
valid only for F v oc y _1 sources and a small colour correc- 
tion has been applied to the "blue" sources, which have sig- 
nificantly different spectra (F v oc y~ 3 ). The reference wave- 
lengths are defined at 6.7 and 14.3 /mi for LW2 and LW3, 
respectively. Conversion from flux density to magnitude is 
defined as m(,j = -2.5 log(f v (6.7/(to)/82.8) and m 1 4 3 A(m 

2.5 log(F v (14.3//m)/18.9), where F v is given in Jy. The ~ 5 % 
responsivity decrease throughout orbit has not been corrected 
for. The ISOCAM astrometry was based on the ISO point- 
ing only. We manually cross-correlated bright sources with 
2MASS positions and found a bulk offset of 7" in DEC and 
1.4" in RA. Correcting for this, the mean deviation between 
ISOCAM and 2MASS positions (for these bright sources) be- 
came 2.5 + 1.4". 

ISOCAM point source extraction and photometry gave a 
total of 186 detections (see TableQ of which 160 have reliable 3 
fluxes in the 6.7 //m band (LW2) and 74 in the 14.3 /mi band 
(LW3). For the 14.3 /mi band we estimate a completeness limit 
of 8 mJy and a limiting sensitivity of 3 mJy, while for the 6.7 
/mi band the survey is complete to 6 mJy and reaches 1 mJy. 
At these levels the sam ple is not co ntaminated by galaxies. See 
the discussion in K aas et al.l ll2004h for details. 

Tabled which is available in electronic form at the CDS, 
lists all the 186 ISOCAM detections with fluxes and uncertain- 
ties for each band (see Sect.[3]for a description of the columns). 
Flux measurements have been obtained in both bands for a to- 
tal of 57 sources. These belong to either of two colour pop- 
ulations, one group of 36 "blue" sources that have colour in- 
dices like normal photospheres, and another group of 21 "red" 
sources that have excess emission at 14.3 /mi (see Fig. |4] and 
Sect. |3] in general for details). Figure \l\ is a sky map of the 
ISOCAM sources. The locations of the IR-excess sources, that 
are also labelled, are shown together with those of the "blue" 
(no IR-excess) sources, as well as sources with either LW2 or 
LW3 fluxes only. The subregions marked in this map show the 
location and extent of various follow-up observations: IRAM 
1.3mm map (large box), Arnica deep Ks imaging (medium 
box), and NOTCam H2 line imaging (small box). The SIRCA 
pointed observations are not indicated in the figure. 



2 CAM Interactive Analysis, a joint development by the ESA 
Astrophysics Division and the ISOCAM Consortium led by the 
ISOCAM PI, C. Cesarsky, Direction des Sciences de la Matiere, 
C.E.A., France. 

3 Detections without reliable flux measurements are generally 
sources observed on the edge of the detector. 



Table 1. The 186 ISOCAM detections in the Ser/G3-G6 com- 
plex. Only available in electronic form at the CDS. 



2.2. 2MASS 

In the region mapped by ISOCAM a total of 1803 sources was 
found in the 2M ASS All-Sky Release Point Source Catalog 
JCutri et al.l20ol . The J (1.24 /mi), H (1.66 /mi), andZT. s (2.16 
/mi) photometry of all these has been used in this paper. 

We cross-correlated the 186 ISOCAM sources with the 
2MASS PSC using a search radius of 7.5" (i.e., 3 times the 
mean positional deviation between ISOCAM and 2MASS for 
bright sources after the global recentering, see Sect. 12. It . Only 
19 ISOCAM detections had no 2MASS counterpart, and of 
these, three are very red and faint in the near-IR, while the re- 
maining had been detected in only one of the two ISOCAM 
bands. In 29 cases we found multiple 2MASS objects within 
the search radius. Disregarding these, the mean positional de- 
viation of all cross-correlations is 3.1 + 1.5" and the maximum 
is 6.5", and this is a measure of the accuracy of the ISOCAM 
positions. None of the 29 ISOCAM sources that are associated 
with multiple 2MASS sources within a 7.5" radius are YSO 
candidates with excess emission at 14.3 /mi, although ISO- 132 
might be a Class II to Class III transition object (cf. Sect.0}. In 
the 28 cases of two 2MASS sources within the search radius, 
we have discarded one of them as "probably not an ISOCAM 
counterpart", based on separation, color, and quality flags. In 
one case of a quadruple (ISO-NH3-1 18), we have discarded 3 
objects. 

2.3. K s band imaging with Arnica at the NOT 

In August 1996, about 4 months after the ISOCAM observa- 
tions, the Arcetri Near-Infrared Camera (Arnica) was used at 
the Nordic Optical Telescope (NOT), La Palma, to obtain deep 
K s band (2.15 /mi) imaging in the most active 8' x 8' subre- 
gion of the ISOCAM field. A 4 by 4 raster mode (256 x 256 
array and 0.54"/pix) was scanned repeatedly (dithering each 
scan), staring at the same position only 1 minute at a time with 
a unit integration time of 1 second, until a total on-source in- 
tegration time of 5 minutes was obtained for the whole map. 
In this way the program images themselves were used to sub- 
tract the thermal sky level. Flat fields were calibrated on the 
decreasing/increasing twilight sky, removing thermal and stray 
light contributions by taking difference images. Standard stars 
taken from the UKIRT Faint standard list (Casali & Hawarden 
1992) and from the Arnica list (Hunt et al. 1998) were observed 
throughout the night. Observ ation and redu ction methods are 
the same as those described in Kaa^ Jl999ah . 

Figure[2]shows the deep Ks band imaging with Arnica of an 
8 x 8 arc minutes sub-region (i.e., 18 %) of the ISOCAM area, 
and contours of the LW2 band (6.7 fim) image are overlaid . 
This area covers the Ser/G3-G6NE ammonia core JClarkll99ll) 
well, and a high extinction region is outlined by the drop in the 
surface density of stars in the Ks band. 

Point sources were detected with daofind (threshold at 5 <x), 
and aperture photometry was made using a small aperture with 
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Fig. 1. Sky map of the ISOCAM sources. Symbol size inversely proportional to LW2 magnitudes. IR-excess sources (filled 
circles), sources without IR-excess or "blue" sources (open circles), only LW2 fluxes (plus sign), and only LW3 fluxes (crosses). 
Coordinates are J2000 epoch. The three subregions indicated are the areas observed with IRAM (large box), Arnica (medium 
box), and NOTCam (small box). The pointed SIRCA observations are not shown. See text for details. 



an aperture correction applied to compensate for flux losses in 
the wings of the stellar profile. A total of 1014 sources (out 
of 1 148 detections) have photometry with ctk s < 0.3 mag in 
the 8' x 8' Arnica map. All are listed in Table |2 (available 
online only). The first column gives the Arnica/NOT identifi- 
cation number from 1 to 1148. The second and third columns 
list the J2000 equatorial coordinates. Columns four and five 
list the Kg magnitude and its error, and the last column gives 
possible other identifications where cross-correlation has been 
made with the ISOCAM map. The coordinates were found by 
registering the mosaic using IRAF tasks ccmap and xy2rd with 
the 2MASS coordinates of 44 cross-correlated ISO and SIRCA 
sources as an input grid to the plate solution. The rms of the 
J2000 coordinates presented in Table[2]are 0.38" in ra and 0.3 1" 
in dec. This relatively large rms is probably because of the op- 
tical distortion in the camera, which was not corrected for. The 
sources with cr Ks < 0.3 mag are in the range from K$ = 7.5 to 
Ks = 17.8 mag, and the survey is estimated to be complete to 
K s a 16 mag. 

Comparing the Ks band magnitudes of Arnica/NOT 
(August 1996) with those of 2MASS (July 2000) for 37 
ISOCAM sources in the range 7.8 < Ks < 14.3 shows overall 
agreement. The difference Ks (1996) - ^(2000) has a median 



Table 2. The 1 148 Arnica K$ band sources detected in a 8' X 
8' subregion of ISOCAM coverage of the Ser/G3-G6 complex. 
Only available in electronic form at the CDS. 



value of -0.005 + 0.07 mag when disregarding six sources that 
have varied by more than 0.2 mag and are therefore variable 
sources. These are: ISO-NH3-84 (a double Class II), ISO-NH3- 
86 (flat-spectrum source), ISO-NH3-94 (Class I or younger - 
the K band emission is probably dominated by the emission 
in the 2.122 ^m S(l) line of H 2 ), ISO-NH3-105 (Class II), 
ISO-NH3-137 (Class II), and ISO-NH3-154 (Class II). Refer to 
Tables 0] and [5] Because the Arnica data was collected only 4 
months after the ISOCAM observations, we use these to calcu- 
late the ff2-7 and o^-u spectral indices. Also, the Arnica dataset 
resolves ISO-NH3-84 into two sources and provides a Ks band 
magnitude for the faint source ISO-NH3-90. 
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Fig. 2. Arnica/NOT /T, band 4x4 mosaic of an 8 x 8 arc minutes part of the ISOCAM area. Contours show the ISOCAM LW2 
(6.7 /im) mosaic and ISO ID numbers are labelled (to the right of the source). 



2.4. V band imaging with SIRCA at the NOT 

The Stockholm Infrared Camera (SIRCA 4 ) was used at the 
NOT, La Palma, 15-19 July 2003. We imaged 12 SIRCA fields 
(56" x 70") within the ISOCAM map in Serp-NH3. This in- 
strument/telescope combination has a sensitivity in the L' band 
(3.8 /im) of 0.2 mJy per square arcsecond, i.e., a limiting mag- 
nitude of L' = 15.2. The detector is a 256 x 320 InSb array, 
and the spatial resolution at the NOT is 0.22"/pix. The images 
were taken with an individual integration time of 0.2 seconds, 
5 reads per chopper position, 5 chopper repetitions, and 20 



4 See lOlofsson Floren. & Lagaed < 19991) and 

Ivan der Bliek et alJ feOOCh . as well as URL 

http://www.astro.su.se/groups/infrared/instrumentation.html for 
more information on SIRCA. 



Table 3. The 46 SIRCA L band detections in 12 pointed obser- 
vations in the Ser/G3-G6 complex. Only available in electronic 
form at the CDS. 



ABBA "noddings" of the telescope, i.e., a total on-source inte- 
gration time of 0.2x5x2x5x20x4 = 800 seconds. The bright 
UKIRT standards H D161903 and HD203856 w ith L' magni- 
tudes 7.01 and 6.84 JCasali & Hawardenlll992l) . respectively, 
were integrated 80 seconds (only two telescope noddings). The 
flatfield was obtained from the target frames by chopping be- 
tween the sky and outside the field, i.e., the cold interior of the 
instrument, which effectively subtracts the dark current. Target 
images are shifted and added for each nodding cycle and flat- 
field corrected. 
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ble and the standards were used instead. The field centred on 
Ser/G3-G6 is shown in Fig.|3] 
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Fig. 3. SIRCA/NOT L' band image of a 70" x 56" area around 
the Ser/G3-G6 group. North is up and East to the left. 



Our pointed observations with SIRCA gave U band pho- 
tometry of 46 sources in 12 fields. All these are listed in 
Table The columns give SIRCA ID number, RA and DEC 
(J2000) positions 5 , L' magnitude and its error, K$ magnitude 
and its error, 2MASS ID, Arnica/NOT K s ID, ISO-NH3 ID, 
and other ID. The K$ band magnitudes are taken from 2MASS 
except for IDs 8, 15, 16, 21, 22, and 34, where they are taken 
from the Arnica/NOT observations listed in Table [2] Positions 
are from 2MASS when a 2MASS source is available, otherwise 
they are ISO positions or positions found from the L' band im- 
age with respect to other sources. For IDs 21 and 22 the posi- 
tions are from the Arnica/NOT images. 

Five repeated measurements of the standard star HD 16 1903 
over 4 nights gave a rms deviation of 0.048 magnitudes. No cor- 
rection was made for airmass, since this would be smaller than 
the variable extinction produced by the presence of dust from 
Sahara at the time of observing. The standards were relatively 
near to the targets, and the instrumental target magnitudes have 
been calibrated with the standard observation nearest in time. 
Magnitudes are calculated with the aperture photometry task 
phot in IRAF, using 20 pix (4.4") aperture radius for standards 
and 5 pix (1.1") for program stars. For the small aperture an 
aperture correction was applied to compensate the loss of flux 
in the wings. The correction was determined individually in 
each image, apart from a few cases where this was not possi- 



5 The RA and DEC (J2000) positions of SIRCA sources are taken 
to be the 2MASS position. If no 2MASS counterpart is found, then 
the position is calculated with respect to other sources in the image. 



2.5. H 2 line imaging with NOTCam at the NOT 

The Nordic Optical Telescope near-IR Camera/spectrograph 
(NOTCam) was used to do narrow-band imaging in the line 
(S(l) 1-0 at 2.122 yum) and a nearby continuum (A c = 2.087 and 
FWHM = 0.02 yum) at the location of the NH 3 core to the north- 
east of the Ser/G3-G6 group. The detector was the 1024 x 1024 
x 18/um Hawaii engineering grade array. Ser/G3-G6NE was 
observed on 28 May 2003. The wide field camera (0.235"/pix) 
was used, and the observations were performed with a ramp- 
sampling readout mode. We integrated for a total of 100s per 
sky position, reading the array every 10s and using the linear 
regression result of these 10 readouts to lower the read noise. 
Altogether 20 dithered sky positions were observed, giving a 
total of 2000 s on-source for each of the two filters. 

Flat fields were observed on the twilight sky creating 
bright -faint pairs for a proper subtraction of the thermal, dark, 
and stray light contributions. Figure ^] shows the H2 line 
(2.122 yum) image of the Ser/G3-G6NE region, and Fig. [H 
shows the narrow-band continuum (A c 2.087 /im) image of the 
same region. The images suffer from some stripy pattern due 
to pick-up noise in the electronics at the time of observing. 
However, it is clear that most of the extended emission is pure 
H2 line emission. Only around the position of ISO-94 and ISO- 
101 (cf. Fig. 1161 faint, extended continuum emission is seen. 

The H2 line image is flux calibrated using t he near-IR stan- 
dards AS 3 3-0 and AS33-1 faunt et alJll998t) . The zeropoint 
difference between the H2 and the K filter was found to be 
2.635 mag. We transform to monochromatic flux through the 
relation F v = 710 x io~°- 4xm (Jy) for a zero'' 1 magnitude star 
at 2.122 yum, assuming the stars' magnitudes are the same at 
2.122 as at 2.2 yum. From the H2 line image of the standards, 
we get a conversion factor of 15.75 yuJy per ADU/s for the im- 
age. The intensities in the knots range from 2 to 46 /Jy/pix, and 
the noise in the background is cr J/tv = 0.5 //Jy/pix. 

2.6. IRAM 30m telescope observations 

A 1.3 mm dust continuum mosaic of part of the ISOCAM 
coverage was taken with the IRAM 30-m telescope equipped 
with the MPIfR 37-channel bolometer array MAMBO-I 
(Kre vsa et alJ ll998) during four nights of observing sessions 
in March 1998. The passband of the MAMBO bolometer ar- 
ray has an equivalent width w 70 GHz and is centred at v e // ~ 
240 GHz. 

The mosaic consists of 18 individual on-the-fly maps that 
were obtained in the dual-beam raster mode with a scanning 
velocity of 8"/sec and a spatial sampling of 4" in elevation. In 
this mode, the telescope continuously scans in azimuth along 
each mapped row while the secondary mirror wobbles in az- 
imuth at a frequency of 2 Hz. A wobbler throw of 45" or 60" 
was used. The typical azimuthal size of individual maps was 4'. 
The size of the main beam was measured to be ~ 11" (HPBW) 
on Uranus and other strong point-like sources such as quasars. 
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The pointing of the telescope was checked every ~ 1 hr using 
the VLA position of the strong, compact Cl ass source F IRS 1 
in the Serpens Cloud Core (good to ~ 0.1". lCuriel et alJl993h : 
it was found to be accurate to better than ~ 3". The zenith atmo- 
spheric optical depth, monitored by 'skydips' every ~ 2 hr, was 
between ~ 0.2 and ~ 0.4. Calibration was achieved through on- 
the-fiy mapping and on-off observation s of the primary calibra- 
tor Uranus (e.g.. lGriffin & Ortonl 19931 and references therein). 
In addition, the Serpens secondary calibrator FIRS1, which has 
a 1.3 mm peak flux density ~ 2.4 Jy in an 11" beam, was ob- 
served before and after each map. The relative calibration was 
found to be good to within ~ 10% by comparing the individual 
coverages of each field, while the overall absolute calibration 
uncertainty is estimated to be ~ 20%. 

The dual-beam maps were reduced and combined with 
the IRAM software for bolom eter- array data ("NIC"; cf. 
Broguier e. Neri & Sieversl l995t). which uses the EKH restora- 
tion algorithm JEmerson""Klein & HaslarrJl9"7^ . 

Figure El shows the IRAM data contour map with four 
bright sources labelled (MMS1, MMS2, MMS3, and MMS4), 
and Table lists the individual sources with fluxes and posi- 
tions. 

2.7. VLA 3.6 cm observations 

Radio continuum observations of three fields covering all four 
millimeter sources (MMS 1 -4) were made at 3 .6 cm (central fre- 
quency 8.46 GHz, total bandwidth 100 MHz) with the NRAO 6 
Very Large Array (VLA) in a mixed B/C configuration on 
22 June 2001. The on-source integration time per field was 
~ 75 min. The amplitude calibrator was 3C286 and the phase 
calibrator was 1801+010. The data were edited and calibrated 
using standard VLA procedures and the corresponding map 
cleaned using the AIPS task IMAGR with natural weighting. 
The synthesised beam size was 2.8" x 2.3" (HPBW) and the 
primary beam (or field of view) of each observation was 5.3' 
(FWHM). The rms noise in each field was ~ 16//Jy/beam, im- 
plying a 5 <x detection threshold of ~ 0.08 mJy/beam in the in- 
ner parts of the fields. A total of eight point-like radio sources 
were detected above the 5 <x level within the primary beams 
of the three observed VLA fields. Table|S]lists their positions 
(accurate to +0.5") and their 3.6 cm flux densities corrected 
for primary beam attenuation. Five VLA sources appear to be 
closely associated with infrared and/or millimeter sources (cf. 
Table [8}. The other sources are probably background extra- 
galactic sources. In addition to these robust VLA detections, 
we have also included a weak, tentative (~ 3.3 cr) radio source 
in Table|8]as it coincides within 1.5" with the 1.3 mm contin- 
uum source MMSl-c. 

3. IR-excess sources 

The results from the different datasets on the ISOCAM se- 
lected sample are summarised in Table ^ which is available 

6 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under a cooperative agreement with the 
National Science Foundation. 



at the CDS in electronic form only. Its first column lists the 
ISO identification number (ISO-nh3-#), then the ISOCAM po- 
sitions (RA and DEC in J2000 epoch) in columns two and 
three. Columns four and five give the Arnica K s band mag- 
nitude and error (K s (96) and <t k ). Columns six to nine list 
the ISOCAM fluxes with uncertainties in mJy (F v (6.1fi), <r 6 7 , 
F v (14.3yu), cr I4 . 3 ). Columns 10 to 13 show the ISOCAM detec- 
tion and photometry flag, the name of the 2MASS counterpart, 
the id of the NOT (Arnica Ks (96)) counterpart, and at last a col- 
umn with possible other identification. We refer to the 2MASS 
catalogue for further information on the 2MASS counterparts. 

In the following subsections we will use various colour- 
magnitude and colour-colour diagrams to extract YSO candi- 
dates based on IR excess. The strategy is as follows: 

- select IR-excess sources from the ISOCAM data 

- compare with the J - H/H - K s diagram from 2MASS and 
search for more candidates 

- combine bands of ISOCAM and 2MASS to search for ad- 
ditional IR-excess sources in the H - Ks /Ks - m-j diagram 

- combine SIRCA L band with 2MASS H and K s to search 
for more candidates 

3.1. The ISOCAM colour-magnitude diagram 

For the ISOCAM sources we present a colour-magnitude di- 
agram in Fig. |4] The colour index [14.3/6.7], defined as 
log(F v (14.3/im)/F v (6.7yum)), is plotted against the flux in the 
14.3 /im filter (magnitude indicated on upper x-axis). This 
colour index can be related to the spectral energy distribution 
(SED)indexo' = -(c/logUF^/^log/i) calculated between 6.7 
and 14.3 jum, which is indicated on the right-hand y-axis. As 
seen in all star formation regions surveyed with ISOCAM, the 
sources tend to separate into two distinct groups: sources for 
which the colour index follows that of normal stellar photo- 
spheres, located around [14.3/6.7] = -0.69 or a = -3, and 
indicated with a dashed line in Fig.|4l and sources with intrin- 
sic IR excess. Among the 57 sources with fluxes in both bands, 
21 have mid-IR excesses. These are all YSO candidates. We 
note that they occupy a broad range in brightness (4-900 mJy 
at 6.7 /mi). In addition, their amount of IR excess at 14.3 /mi is 
apparently independent of brightness. 

There are 16 sources detected at 14.3 /mi only and plotted 
as upper limits (small crosses) in Fig.|4] Only 8 of these have 
2MASS counterparts, all with good photometric quality in K. 
The SED index taken between 2 and 14 //m (cf. Sect.0Ji gives 
-1.01 < cfj R 14 < 1-2, values corresponding to an amount of 
IR excess typical of Class I and II types of YSOs. These are: 
ISO-60, 136, 143, 153, 166, 175, 181, and 182, of which ISO- 
143 is extremely red and outside of the plotted region in Fig.|4] 
We add these 8 to our list of YSO candidates. The remaining 8 
might be spurious sources. 

There are as many as 103 sources detected at 6.7 /mi only. 
For 97 there is a 2MASS counterpart, and all but 6 of these have 
good photometric quality in all three bands. These sources are 
discussed in Sect. 13. 31 
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Fig.4. The colour index [14.3/6.7], defined as 
log((Fy(14.3//)/Fy(6.7jU)), is shown on the y-axis and the 
flux at 14.3/im (in Jy) is shown on the x-axis. Sources without 
IR excess tend to line up around [14.3/6.7] = -0.66 or a = -3 
(dashed line). The 21 sources with mid-IR excesses (circles) 
are spread out over a large range in brightness. These are all 
YSOs. Lower limits on colour are given for sources detected at 
14.3 fim only (small crosses). 



Fig. 6. The H-K s /K s -m 7 diagram for 151 ISOCAM sources. 
The reddening vector indicated for an AO star (thin continuous 
line) is found by a linear fit to the blue ISOCAM sources (open 
triangles). The dotted lines indicate the reddening band of gi- 
ants, supergiants, and main-sequence stars. The red ISOCAM 
sources (open circles) are practically all to the right of the red- 
dening band. The 12 additional sources with IR excess found 
from this plot are marked as large asterisks. 



Aj. ~ A la 




□ 



Fig. 5. The J-H/H-K s diagram for the ISOCAM sources with 
2MASS counterparts (dots). The red ISOCAM sources (open 
circles) merge together with the blue ISOCAM sources (open 
triangles) along the reddening band of normal stars. Seven 
ISOCAM sources appear to have excess emission in the H-Ks 
index (open squares), but they are all flagged 2MASS sources. 
Twenty faint 2MASS sources (crosses) with apparent IR ex- 
cesses are discussed in the text. 

3.2. The 2MASS J - H/H - K s diagram 

Figure0shows the / - H/H - K s diagram for all the ISOCAM 
sources with 2MASS counterparts (dots). The blue ISOCAM 
sources are marked with open triangles. The red ISOCAM 
sources found to be YSO candidates in the previous section 
are marked with ope n circles. The loci of giant, supergiant, and 
main-sequence stars (Koornneef 1983) are indicated with bold 
curves. We have calculated the slope outlined by the ISOCAM 
blue sources to be 1.82±0.05. This should be a good indica- 
tor of the cloud extinction in this diagram, and the redden- 



ing vector for an AO star is shown by the dashed line. For the 
2MASS JHK S filter pass bands at AA 1.235, 1.662, and 2.159, 
respectively, this slope is in accordance with the A] <x A~ 
parametrization of the NIR extinction law ( Whittet 1988). This 
same opacity index /3 — 1 .7 ±0.4 was found for Se rpens in the 
JtoH wavelength region bv lFroebrich et alJd2004l) . also based 
on 2MASS data. 

Figure|5]shows that the majority of the sources are clustered 
along the reddening band. To distinguish IR excess sources 
from reddened sources in the J — HI H-Ks diagram, they must 
be located to the right of the reddening line by at least 1 <x of 
the individual uncertainties in the H-Ks index. Only 10 of the 
29 mid-IR excess sources have a detectable IR excess in the 
near-IR, as measured from the J - H and H - Ks colours. Thi s 
is in line with the findings of, e.g.. lKaas & Bontempsl (EoOO). 
that only half or even less of the IR-excess population found 
from mid-IR colours is sampled by using near-IR colours. 

There are 7 ISOCAM sources, for which no IR excess could 
be found from the ISO observations, which seem to have ex- 
cesses in their H-Ks color (see Fig. [5}- We note that all of 
these are upper limits in at least one of the 2MASS bands and 
should not be considered YSO candidates. Two of them are 
blue ISOCAM sources (ISO-59 and 1 16), one is independently 
found to have IR excess (ISO-64, see next section), but the na- 
ture of the remaining 4 cannot be assessed (ISO-31, 53, 91, 
128). 

In the ISOCAM field there are 1221 2MASS sources with 
photometric quality flag 7 A to D in all three JHKs bands. 
Among these there are many faint sources not detected by 
ISOCAM with apparent near-IR excesses. Because these are 
faint sources at the survey limit, we consider those located to 



7 See the 2MASS catalogue description at 
http://www.ipac.caltech.edu/2mass/releases/allsky/doc/ for details. 
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the right of the reddening band by at least 2 cr of the individ- 
ual uncertainties in the H - Ks index. The resulting 20 sources 
(see Fig. [5} occupy a range of J - H and H - K s colours typi- 
cal of pre-main sequence stars. As suggested by the referee, we 
cross-correlated these with the on-line c2d Spitzer photometry 
table for Serpens, and found IRAC counterparts for 19 sources. 
None of these have IR excess in any of the Spitzer IRAC bands, 
however, and therefore they cannot be proposed as bona-fide IR 
excess YSO candidates. The apparent excesses in H - Ks are 
not understood and may mostly be the result of statistical fluc- 
tuations in the 2MASS data, since none of the 20 sources is 
located to the right of the reddening band by more than 5 cr of 
the individual errors in the H - Ks index. More data is needed 
to establish the near-IR colours of these sources. 



3.3. The H - K s /K s - ra 7 diagram 

As seen in Fig. [5] the J - H/H - Ks diagram does not extract 
nearly half of the IR-excess sources identified with ISOCAM. 
The IR excess does not separate well from the reddening in 
the J - H and H - K s colours, and the blue and red ISOCAM 
sources tend to merge together al ong the reddening band. As 
shown in, e.g.. iKaas et alJ 12004), the H - Ks /Ks - m-j dia- 
gram better distinguishes IR excess from reddening. By using 
this diagram we can study the 103 ISOCAM sources for which 
there are 6.7 //m fluxes only (see Fig. [6}. The reddening vector 
indicated is empirically found from this dataset. The slope 1 .09 
+ 0.01 is found by a linear fit to the blue ISOCAM sources. 
For c omparison, the sl ope was 1.23 for the Serpens Cloud 
Core JKaas et al.lEo 04). The width of the reddening band of 
giants, super giants and main-sequence stars, whose intrinsic 
colours are outlined by the bold face curves, is indicated. The 
red ISOCAM sources are all, except two (132 and 167), found 
to the right of the reddening band. We find a total of 26 IR- 
excess objects from the H - Ks /Ks - mj diagram, i.e., those 
who are located to the right of the reddening band of normal 
stars by more than 1 cr in the Kg - mj colour index. Twelve 
of these are new, i.e., not selected in the colour diagrams dis- 
cussed above. 



3.4. The H-K s /K s -L diagram 

The H-Ks /Ks - L' diagram is shown in Fig.0 We plot the 39 
sources available in these colours, and overplot red ISOCAM 
sources, blue ISOCAM sources, and additional IR-excess 
sources found from the H - Ks /Ks - mj diagr am. The intrin- 
sic co lours of main sequence and giant stars ( Bessell & Brett 
1988) are indicated, too. The redden ing vector of a n AO star is 
drawn applying the Aj oc law dWhittedll988h again, giv- 
ing a slope of 1 . 1 8. The reddening band for normal stars is also 
indicated, and IR excess sources must be located to the right 
of this line by more than one cr of the individual errors in the 
Ks - L' colours. We note that the distinction between the red 
and blue ISOCAM sources is not so clear in this diagram as in 
the H - Ks /Ks - ni-j diagram, but significantly clearer than in 
the J-H/H—Ks diagram. There are no blue ISOCAM sources 
in the IR-excess locus, but there are five red ISOCAM sources 
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Fig. 7. The H - K s /K s - L' diagram for 39 sources (dots). We 
show by overplotting symbols the red ISOCAM sources (open 
circles), blue ISOCAM sources (open triangles), and the ad- 
ditional IR-excess sources identified in the H - Ks /Ks - mj 
diagram (large asterisks). The reddening band of normal stars 
is indicated (dotted line) as is the reddening line for an AO star 
(thin continuous line). Four new IR excess sources are identi- 
fied. 

in the locus of reddened normal stars. In total 1 1 sources have 
IR excess in the Ks - L' colour, of which four were not previ- 
ously identified as IR-excess objects. None of these four is in 
the ISOCAM selected sample of Table but two of them are 
the objects G4 and G5, which are very close to G3 (ISO-82) 
and G6 (ISO-78), respectively. The third object is in the same 
area just SE of G4, and the fourth object is a faint star in the 
SIRCA field of ISO-167. 

3.5. Extracted YSO candidates 

From the 186 ISOCAM detections we have extracted 29 YSO 
candidates from the ISOCAM data, and 12 additional from the 
combination of ISOCAM and 2MASS using the H - K s /K s - 
m-j diagram. This gives a total sample of 41 IR-excess YSOs 
within the ISOCAM selected sample. In addition, we find 4 
IR-excess objects from the H -Ks /Ks —LI diagram. Summing 
up, this gives a total of 45 YSOs. 



4. The YSO population 

According to the current empirical picture of young stel- 
lar evolution, low mass YSOs go through evolution- 
ary phases with different Spectral Energy Distribution 
(SED) characteristics: Class 0, Class I, (flat-spectrum,) 



1987 



1993 



Class II , and C lass III sources dAdams. Lada & Shul 
[ Ladal Il987t L \ndre. Ward-Thompsom&Barsoiivr 
Andre & Montmerle 1994). 

The Class sources are in the deeply embedded main 
accretion phase and have circumstellar envelope masses 
larger than their central stellar mass. The defining char- 
acteristics of a Class is: strong, centrally condensed 
dust continuum emission in the submm, very little emis- 
sion shortward of 10 fim, and powerful jet-like outflows 
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Fig. 8. The spectral index a 2 ~ 7 versus ff^ 14 or equivalently, 
the K - m-i/K - mi4 colour diagram. The locations of Class I, 
Class II, and Class III sources and field stars are indicated, as 
well as the flat-spectrum source locus. 



jAndre. Ward-Thompson. & Barsonvll2000l) . Thus, we do not 
expect to sample these sources in the mid-IR, although 
ISOCAM did detect a few Class sources, and Spitzer with its 
improved sensitivity has detected many more. The dust contin- 
uum mapping at 1.3 mm with IRAM should sample the Class 
sources, and we refer to Sect. [6] for a discussion of these re- 
sults. The Class I, fiat-spectrum, and Class II sources, however, 
have IR SEDs that are reasonably well sampled by flux mea- 
surements at 2.2, 6.7, and 14.3 fim. 

4.1. Classifying the IR-excess sources 

We have calculated the SED index between 2.2 and 14.3 ^m, 
q?p 14 , whi ch is close to the index used by lLada & wTl king 
(1984 and IT ^daN 19871) to define the three classes, I, II, and 
III, originally . We use the updated IR spectral classific ation 
scheme ( Andre & Mon tmerld 1 1994 iGreene et all 1 1994 and 
define YSOs with a 2 IR u > 0.3 as Class I sources, YSOs with 
a 2 R 14 in the interval [-0.3,0.3] as flat-spectrum sources, and 
YSOs with a 2 R 14 in the interval [—1.6, -0.3] as Class II sources. 
IR-excess sources with a] R lA < —1.6 are tentatively defined as 
Class III sources. We note that because our YSO sample is se- 
lected from IR excesses most Class III objects are missed. Only 
a few with a small detectable IR excess are captured, probably 
those objects in transition between clasees II and III. Class III 
sources with no IR excess cannot be distinguished from field 
stars using our data. Additional tools such as X-ray surveys 
must be applied to sample the Class III population. 

Figure [8] shows the a 2 R 7 index versus the a 2 R 14 index with 
the corresponding K - mq and K - mu colours in magni- 
tudes given on the opposite axes. The loci of the different YSO 
classes are indicated. Sources with no IR excess line up in the 
lower left of the diagram along the extinction vector. The vector 
is taken from lKaas et al. (2004) and has a length corresponding 
to ~ 60 mag of visual extinction. Most sources are found with 
A v about 10-20 mags, while one object (ISO-88 = HD 170377) 
has no extinction and is perhaps a foreground star or a Class III 



cluster member at the front side of the cloud. According to 
SIMBAD it is a bright V = 8.96 mag optical star with spectral 
type K2. Sources with IR excess mainly occupy the Class II lo- 
cus, while there are also a few in the Class I, flat-spectrum, and 
Class III loci. Compared to the Serpens Cloud Core JKaas et alJ 
2004) this region has fewer Class I and flat-spectrum objects. 
The reddest Class I, however, has a 2 R 14 = 2.94, which is redder 
than any of the Class Is in the Serpens Cloud Core. Both very 
red objects (ISO-90 and ISO-94) are located in a core of active 
star formation near the dense mm core MMS2 (see Sects. I6~T1 
and 16. 2> . At least for one of these (ISO-94) the emission in the 
K-band is dominated by line emission. 

We note the almost linear relationship between the two 
SED indices a 2 R 7 and a] R u in Fig. HO As in lKaas etal] J2004) 
we use this relationship to classify objects from the a 2 R 7 in- 
dex only, allowing us to classify objects not detected at 14.3 
jt/m. From Fig.[S]we see that the colour [14.3/6.7] is essential 
to distinguish Class III objects with IR-excess from reddened 
stars. It is also evident that for very high extinction there is a 
risk of confusing Class II objects with reddened stars unless the 
[14.3/6.7] colour is known. Nevertheless, with this caveat in 
mind, we classify 12 more sources than those plotted in Fig. [8] 
Thus in the sample of IR-excess YSOs, we define as Class Is 
those that have a 2 R > 1.2 (K - m 7 > 4.8), Class lis those that 
have a 2 R 7 < -0.25 (K - m-; < 3.2), and flat-spectrum sources 
as those in between. 

There are, however, 6 IR-excess sources with neither of 



the two SED indices (a 2 R 14 



and a 2 IR ) available. In the follow- 



ing we will classify these sources on an ad hoc basis. Four of 
them are SIRCA sources (id 9, 1 1, 13, and 44 in TableGJ with 
no ISOCAM detections. Plotting the various YSO classes in 
the H - K/K - L' diagram we see that these four are located 
in the region occupied by the Class II sources in our sample. 
Therefore we tentatively classify them as Class lis. ISO-62 is 
a red ISOCAM source without a 2MASS counterpart located 
in the dense Ser/G3-G6SW ammonia core. Judging from the 
2MASS sensitivity we estimate that ISO-62 has a 2 IR u > 2 and 
we tentatively classify it as a Class I object. ISO- 101, located 
in the dense Ser/G3-G6NE ammonia core, is a red ISOCAM 
source not detected in the near-IR continuum. The NOTCam 
images show extended emission in the H2 line at 2.122 fim and 
no emission in the continuum filter at 2.087 //m. Its rather blue 
[14.3/6.7] colour probably also reflects a H2 line emission (cf. 
Fig.ll3i. Some faint extended emission in the L' band is barely 
detected at this position (see Table[3J. The lower limit estimate 
of a 2 R 14 first suggested that ISO-101 could be a Class I or flat- 
spectrum source, but considering the above discussion we do 
not count this object as a continuum source (see discussion in 
Sects.lOandlO. 

In total, the IR-excess population we have found in the 
Ser/G3-G6 complex consists of 5 Class I sources, 5 flat- 
spectrum sources, and 31 Class II sources. In addition, we see 
3 Class Ills with some small IR-excesses, noting that this class 
is not sampled by our IR-excess selection criterion. 

The Class I/Class II number ratio of 5/3 1 is more or less 
typical for sta r-forming regions a nd compares with 16/123 in 
p Ophiuchi (Bontemps et al. 2001) and 5/42 in ChamaeleonI 
(Kaas 1999b). This is in strong contrast to the large fraction 
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Table 4. The 5 Class I and the 5 flat-spectrum so urces. Cross- 
corre lated with Spitzer YSO candidates from lHarvev et alJ 
(2006). 



ISO 


Other ID 


2-1 
U IR 


a- 14 


AK 


90 


TRAS 18?fS5+00?8 Snit7er-1(S 


4.42 


2.94 




O 1 1 


opllZci - 1 / 


J. JO 


Z.jo 


U.Of 


62 2 




>3 


>2 


n.a. 


143 3 






1.20 


n.a. 


64 


Spitzer-5 


2.93 




n.a. 


86 




0.31 


0.27 


-0.74 


60 


Spitzer-3 




0.02 


n.a. 


182 






-0.04 


n.a. 


153 






-0.14 


n.a. 


121 




1.06 




n.a. 



1 H 2 line emission probably dominates the K band flux. 

2 Lower limits on ff 2 ^ 14 , but a^ 14 = 1.54 . 

3 Not located within the 1.3 mm emission, and its Class I nature is 
questionable. 



of Cla ss I sources found in the Serpens Cloud Core jKaas et all 
2004). 

4.2. Class I and flat-spectrum candidates 

There is a total of 10 IR-excess protostar candidates listed in 
Table |4] 5 Class Is and 5 flat-spectrum sources. We note that 
the classification is tentative in the sense that only lower limits 
on the SED indices can be given in the case of ISO-62. Also, 
the indices of ISO-94 are lower limits since the K-band flux 
is dominated by H2 line emission. These sources may be even 
younger than we have indicated. Nevertheless, in the following 
we will use our tentative classification. 

Two of the Class Is are located in the ammonia core Ser/G3- 
G6NE (ISO-90 and 94), and another two (ISO-62 and 64) in the 
core Ser/G3-G6SW. As shown in the IRAM 1.3 mm continuum 
map in Fig El all but one Class Is are located in the high den- 
sity regions outlined by the 1.3 mm continuum emission. The 
exception, ISO- 143, lies to the NE of the two ammonia cores 
in the vicinity of ISO-141, It is extended in both ISO bands and 
not detected as a point source at 6.7 /mi (therefore not visible 
in Fig. |2jl. Its Class I nature might be questionable. The flat- 
spectrum sources are partly located inside the cores and partly 
outside. 

4.3. The Class II and Class III sources 

Our sample of pre-main sequence stars in the Ser/G3-G6 com- 
plex comprises the Class II objects listed in Tableland the few 
Class Ills our selection method is capable of seeing, listed in 
Table [6] As emphasized in Sect. 14. II our YSO selection crite- 
rion is IR excess, and therefore we do not sample the Class III 
population except for occasional transition objects between 
Class II and HL 

As shown in Fig.[Q]the spatial distribution of the Class II 
sources is relatively scattered, apart from the small cluster of 
Class lis at the location of CoKu-Ser/G3-G6. We note that 



Table 5. The 3 1 Clas s II sources. The cro ss-correlated Spitzer 
YSO candidates from lHarvev et al.l (|2006) are indicated in the 
last column. 



ISO Mj Aj AK Other ID 



7 


4.72 


3.51 


n.a. 




26 


5.67 


4.44 


n.a. 




74 


3.50 


4.88 


0.18 




78 


3.06 


1.74 


0.05 


CoKu-Ser-G6, NOT- 158 


80 


6.06 


2.58 


n.a. 


SIRCA-5,Spitzer-10 


81 


2.80 


6.28 


0.04 


SIRC A-7,NOT- 1 62,Spitzer- 1 1 


82 


2.31 


1.63 


-0.11 


CoKu Ser G3, NOT-201 


84' 


5.03 


3.74 


-0.72 


SIRCA-15 & 16, Spitzer-14 


96 


4.76 


3.68 


n.a. 


SIRCA-23 


100 


3.90 


5.73 


n.a. 


SIRCA-28 


102 


5.67 


3.64 


n.a. 




103 


-0.76 


6.70 


n.a. 


SIRCA-30 


105 


0.55 


3.47 


0.50 




122 


3.41 


4.03 


n.a. 




136 


7.00 


4.39 


n.a. 




137 


-0.28 


4.51 


0.49 


Spitzer-21 


139 


2.42 


5.41 


n.a. 




140 


4.28 


4.49 


n.a. 




141 


4.43 


0.22 


-0.07 


SIRCA-41,NOT-898 


147 


3.80 


4.55 


n.a. 




154 


5.70 


3.72 


0.25 


Spitzer-23 


156 


4.78 


3.94 


n.a. 




166 


6.53 


3.48 


n.a. 


Spitzer-24 


175 


9.54 


0.19 


n.a. 




176 


4.46 


0.78 


n.a. 




178 


3.40 


5.18 


n.a. 




181 


6.32 


2.49 


n.a. 






4.26 


2.04 


-0.28 


CoKuSerG5, Spitzer- 12 




3.23 


1.18 


-0.10 


CoKu Ser G4 




7.10 


0.79 


-0.14 


SIRCA-13 




5.40 


3.97 


n.a. 


SIRCA-44 



1 Double source resolved by SIRCA and NOT/Arnica (NOT-239 & 
247), but not resolved by ISO and 2MASS. 



all four objects CoKu-Ser/G3-G6 , earlier found to be T Tauri 
sources with strong Ha emission JCohen & Kuhill97^) are in- 
dependently classified here as Class II sources. For the Class II 
sample we present Mj and Ay as calculated in Sect. 15.11 
in Table [5] The fourth column represents the variation in 
the K s magnitude between the 2MASS data (Jul-2000) and 
the Arnica/NOT (Aug- 1996) dataset in the sense ^(1996) - 
Ks (2000). Among the 12 Class lis where Ks band photometry 
is available for both epochs, as many as 6 (or 50 %) have varied 
by more than AK = 0.2 mag. 

4.4. Comparison with recent Spitzer results 

At the finalization of this paper we be came aware of the re- 
cent publication of Harv ey et alJ (2006) who list Spitzer mid- 
IR fluxes for 24 YSO candidates found in a 12' x 12' area cen- 
tred on the Ser/G3-G6 complex, which they refer to as Serpens 
Cluster B. We have cross-correlated these with our results and 
find that we have 4 Class I and 7 Class II sources in common, 
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Table 6. The 3 Class III sources with some IR-excess. Note that 
our study does not sample the Class III sources. 



ISO 


2-7 


a 2 -' 4 


J-H 


H-K 


K-L' 


AK 


2 


-2.06 


-1.86 


1.619 


0.638 


n.a. 


n.a. 


132 


-2.36 


-2.10 


1.511 


0.735 


n.a. 


n.a. 


167 


-2.19 


-1.86 


1.628 


0.691 


0.435 


n.a. 



while lHarvev et al 1 d2006l> classify one blue ISOCAM source 
as a Class II and find 1 1 YSO candidates that go undetected by 
ISOCAM. The Spitzer source number is added under Other ID 
in both Tables @] and [5] Comparing the SED indices with the 
Spitzer results, the classification is in agreement for all but one 
source, ISO-64 (Spitzer-5), where the Spitzer data suggests it 
is a Class II. The ISOCAM 6.7 /im flux is about twice as high 
as the interpolated value between the Spitzer bands at 5.8 and 
8.0 fim. This could be due to strong H2 line emission in the 6.7 
fim band or also to source variability. 

We also note that both ISO-62 and ISO- 101 are clearly de- 
tected in the IRAC images, but not l isted among the Spitzer 
YSO candidates of lHarvev et al 1 ll2006l) . probably because of 
their extended nature. 

5. Luminosity and mass function 

5.1. Stellar luminosities 

We do not attempt to estimate luminosities and masses for the 
protostar candidates (flat-spectrum and Class I sources). In this 
section we concentrate on the Class II sources. It is probably 
justified to assume that our Class II sample in Table [5]is dom- 
inated by Classical T Tauri Stars (CTTS), since there are only 
very few examples of weak-lined T Tauri Stars (WTTS) with 
the amount of IR excess compatible with a Class II designation. 

To estimate the reddening of each object we assume 
an intrinsic colour of (/ - H)n _ 0.85, i.e., the me- 
dian value found for CTTS JStrom. Strom & Merrilll 1 19931 
iMever. Calvet. & Hillenbrandll 19971) . The extinction in the /- 
band is then found from the relation Ay = 2.53 x [(/ - H) - 
(/ - H)o], obtained from the oc /T 1 7 extinction law (cf. 
Sect. 13. 21 . The absolute /-band magnitudes are My = mj-Aj- 
DM, wher e the distance modulus is DM = 6.76 f or the distance 
of 225 pc JStraizvs. Cernis. & Bartasiute II2003T) . Table lists 
the Class II YSOs with extinction and absolute /-band magni- 
tudes. 

Apart from the photometric uncertainty in /, which is 
+0.05 mag on the average, the uncertainties introduced in My 
arise from 1) the spread in the intrinsic (/ - H)o colour of 
our targets, 2) the uncertainty in the distance estimate, and 
3) the amount of veiling in the /-band. The spread around 
the median value of (/ - H)o for CTTS is of the order of 
+0.15, which translates directly into +0.38 mag in the esfi- 
mat e of Ay. The quoted uncertainty in the distance of +55 
pc JStraizvs. Cernis. & B artasiute1 l2003l) gives an error in the 
distance modulus of magnitudes. The average amount of 
/-band veiling for 22 C TTS in Taurus-Auriga was found by 
Folha & Emerson (1999) to be as high as r = 0.57. About 




M, 



Fig. 9. The My luminosity function for the Class II sources The 
range of young BDs and young stars is indicated by the sub- 
stellar limits (vertical lines) for the two ages 1 and 10 Myr. 




A v = 10 mag 



Fig. 10. The My vs. / - Ks colour magnitude diagram for the 
Class II sources. The 2 Myr isochrone of Baraf fe et alJ ([1998) 
and reddening vectors of length Ay = 10 mag are drawn. Note 
that the J— Kg colour is dereddened according to the deredden- 
ing scheme explained in the text. 



44 % of the CTTS have this or a higher amount of veiling. 
This means that in nearly half of the cases, as much as 36% or 
more of the flux in the /-band (r/(l + r)) is non-photospheric. 
The errors introduced by not correcting for the veiling can thus 
be of the order of 0.4 mag in the direction of overestimating 
the source brightness. On the assumption that these sources of 
errors are independent, we have a total uncertainty (added in 
quadrature) in the estimate of My of ^ ^ magnitudes. 

The My luminosity function of the pre-main sequence sam- 
ple is shown in Fig. [9] A binsize of 0.8 mag was chosen to 
account for the possible uncertainties in the My estimates dis- 
cussed above. The vertical lines indicate the substellar limit for 
the two extreme assumptions on age: 1 Myr and 10 Myr 8 . 

The My distribution is relatively flat over the whole sam- 
pled range. The completeness limit of the Class II sample is 



8 The age of CTTS is currently believed to be of the order of a few 
Myr, but t here is some eviden ce of longlived ~ 10 Myr old circumstel- 
lar disks lLawson et all2002LlLvo et all2003l) . 
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based on the completeness limit of the mid-IR survey, which 
is 6 mJy for the 6.7 /mi band (cf. Sect. 12. It . We have made a 
linear fit between the 6.7 fim flux and My for the 23 Class II 
sources that have both values. The 6 mJy completeness limit 
was found to correspond to My = 5.6 + 0.2 mag. 

In Fig. EH we plot My versus J - K s for the 31 Class II 
sources. The 2 Myr isochrone from the evolutionary models of 
Bar affe et alJ Jl998l) is drawn with the reddening vectors indi- 
cated for the lowest and highest mass included in the model, as 
well as the substellar mass limit at 0.075 M Q . 

The J - Ks colour has been corrected for reddening as dis- 
cussed above. Because of intrinsic IR-excess and possibly also 
insufficient dereddening in some cases, practically all sources 
fall to the right of the 2 Myr isochrone. The sample comprises 
objects in the mass range from above 1.2 M to 0.07 M Q (one 
object reaching to 0.02 M ), assuming an age of 2 Myr. 

5.2. Candidate young brown dwarfs with disks 

The border between stars and brown dwarfs (at ~ O.O75M ) 
is indicated both in Figs. [9] and [K)] We note that between 1 
and 5 of the 3 1 Class II sources in Table [5] have substellar 
masses, judging from Fig.[9]and the age range 1-10 Myr. This 
same fraction of BDs among ISOCAM selected Class II ob- 
jects (i.e., ~ 20 % for a ges of a few Myr) was also found in th e 
Chamaeleon I cloud (Olof sson et alJll998t Persi et alJ EoOO). 
the p Qphiu chi cloud feontemps et al.ll2001 ), and the Serpens 
Cloud Core JKaas et al.l2004 . 

All the sources in our YSO samples have IR-excess, of 
which 31 sources are designated Class II objects. Class lis have 
a mid-IR SED compatible with circumstellar dust distributed 
in a disk. The substellar objects in our study are therefore most 
probably young BDs with circumstellar disks. We note that 
since the selection criteria in our study is IR excess, we do not 
sample any BDs without disks. 



5.3. The Serpens Class II mass function 

The 3 1 Class II sources from this paper are added to the sample 
of 43 Clas s II sou rces in the Serpens Cloud Core region from 
iKaas et a"fl J2004I) to improve the small number of statistics for 
the Serpens Class II luminosity function (LF). The two samples 
are from two different embedded young clusters separated by 
a projected distance of 45 arc minutes (~ 3 pc), but probably 
located in the same cloud system and at the same dista nce (225 
pc adopted from lstraizvs. Cernis. & Bartasiute 12 003). 

For the 3 1 Class II sources of this paper (cf. Table |5Jl we 
have calculated the stellar luminosity L* from My using t he re- 
lation log L+ - 1.49 - 0.466 * My found by iBontemps et all 
il200ll) . The total Serpens sample now consists of 74 Class II 
sources. Figure ^2 shows the observed luminosity function 
(LF) of the combined sample of Class II sources together with 
model LFs calculated for co-eval formation for a set of four 
ages: 0.5, 1, 2, and 5 Myr. 

We assume as a first approximation co-eval formation of 
the Class II sources. Though admittedly a simplification, this 
assumption may be viable as we are treating a subset of cluster 





1 0.5 Myr ' 









Fig. 11. The histogram shows the observed LF of Class II 
sources from the Ser/G3-G6 region (this paper) and the Serpens 
Cloud Core JKaas et al.l2004 as the number of sources per bin 
vs. log L. Overplotted are the synthetic LFs of a coeval popu- 
la tion calculated on the basis of the pre-m ain sequence tracks 
of lD'Antona & Mazzitellil ( ll998l) . Both the lScalol ill 9981) three- 
segment power-law IMF (solid line) and the iRroupal d200ll) 
three-segment power-law IMF (dashed line) trace well the ob- 
served LF for a coeval population of ~ 2 M yr down to the com- 
pleteness limit (vertical dotted line). The Salpeterl (ll955!!) IMF 
(dashed-dotted line) has been included for reference. 



members constrained evolutionarily: Class II sources are cur- 
rently believed to have ages of a few Myrs. Also, it can be ar- 
gued that star formation ma y proceed in seve ral bursts giving 
rise to generations of YSOs (Kaa s et al.l2004l) . 

For easy comparison we sh ow the model LFs here cal- 
culated as in IKaas et all (|2004l) using the pre-ma i n sequ ence 
evolutionary models of iD'Antona & Mazzitellil dl998l) and 
applying three different IMFs. Here we show the Salpeter 
IMF (Salpeter 1955), the Scalo three-segment power-law IMF 
([ Scaloll 19981). and the Kroupa three-segment power-law IMF 
kroupdl200il). W ith the IMF of the form dN oc m r J(logm), 
fhe lSalpeteiHl955l) IMF 9 has T = -1 .35. 

In Fig. ^2 the LFs produced by these IMFs for the four 
ages mentioned above are shown superposed on the observed 
LF, i.e., as the number of sources per bin versus logL. The 
completeness limit at L ~ O.O8L is also given (cf. Sect. 15- 1 1 . 
The normalization has been made to the number of sources 
above the completeness limit. The bin width in the histogram 
is c/ log L = 0.325, based on a factor of two uncertainty in the 
luminosity estimate. 

The incr eased sample in this work confirms the findings in 
dKaas et al J 120 04) that the Serpens Class II LF is compatible 
with coeval formation ~ 2 Myr ago and an IMF of the three- 



r = 

I Mr. 



The lScalol < 19981) IMF has T = -1.3 for 10M o < m < 100M o , 
= -1.7 for 1Mb < m < 10M o , and T = -0.2 for 0.1M o < m < 
The lKroupalfcOOll) IMF has T = -1.3 for m > 0.5M o , T = 



-0.3 for 0.08M G 



sec 



< m < 0.5M H , an d f = 0.7 for m < .08M, „ 
B onnefL^areori^^ZmneckeJ feood) . In lKaas et alJl2004l) we instead 
used KrouDa, Tout & Gilmore 1 1993). 
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Fig. 12. The IRAM 1.3 mm continuum map with the location of the YSOs classified in this paper marked. The contour levels 
go from 30 to 450 mJy/1 l"-beam in steps of 60 mJy/1 l"-beam. The rms noise is 9 mJy/ll"-beam in the best parts of the image 
and typically around 15 mJy/1 l"-beam on average. Class I sources (red circles), fiat-spectrum sources (red crosses), and Class II 
sources (blue filled squares) are indicated. Th e four bright mm sources are labelled MMS1, MMS2, MMS3, and MMS4. The 
0. 12 pc clustering scale for protostars found in lKaas et alJ J2004i) is shown. 



segment power-law type all the way down to the completeness 
limit at L ~ O.O8L , which corresponds to M ~ 0.1 5M Q for 2 
Myr. In the plots we have extrapolated the Scalo IMF and the 
Salpeter IMF to masses below which these IMFs are defined, 
i.e., beyond O.1M and O.4M , respectively. We note that even 
the Salpeter IMF explains the data at high luminosities. 

The peak in the observed LF at L ~ O.O9L coincides 
with the completeness limit and cannot be used to constrain 
the age. An age dependent peak in the LF is expected as a 
result of the piling up of stars in a given luminosity bin be- 
cause of the slowing-d own effect deuterium burning has o n 
YSO contraction (Zinn ecker. McCauehrean & Wilkind fT993). 
Nevertheless, Fig.^Jshows that neither of the three IMFs pro- 
duce compatible LFs for ages much younger or older than 2 
Myr. The model LFs for coeval populations that formed 0.5 
Myr and 1 Myr ago have more sources than observed in the 
luminosity bins around 3L and O.3L , respectively. For ages 
above 2 Myr the model LFs produce fewer sources than ob- 
served for luminosities above O.7L . This tendency starts al- 
ready at an age of 3 Myr, but for 5 Myr the deviation is larger 



than the error bars for several consecutive bins in the observed 
LF. 

We thus conclude that an age of 2 Myr is a good fit for the 
Class II population, although co-eval star formation in a strict 
sense is unlikely. We cannot exclude models where star forma- 
tion proceeds in a more continuous way. However, we find that 
continuous SF over intervals shorter than 2 Myr is incompatible 
with the observed LF. A single burst of star formation that took 
place 2 Myr ago can fully explain the data. Bursts of short du- 
ration produce model LFs that resemble the co-eval LFs, while 
increasing the burst duration smooths out the structure in the 
LFs. 

To distinguish between the IScalol and the Kroupa 

J2.00 lh IMFs for the Serpens Class II sources we need to extend 
the completeness limit from O.O8L to O.O4L . This should be 
feasible with the deeper mid-IR surveys from Spitze r. For the 
more n earby p Ophiuchi region (d = 140 nc). iBontemps et all 
l l200ll) found the IMF of Class II sources (complete to 0.03L G ) 
to be well fitted by a two-step power-law with indices T = 
-1.7 for m > O.55M and F = -0.15 for the whole in- 
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Fig. 13. NOTCam H 2 line image (2.122 jum) of Ser/G3-G6NE. 
Pure line emission features are arrowed (cf. the continuum im- 
age in Fig. I14> . The locations of the two bright mm sources 
MMS2 and MMS3 are marked with plus signs. The Class I can- 
didates ISO-90 and 94, as well as ISO- 101, probably a hot spot 
in the jetlike emission, are encircled. 

terval 0.055M Q < m < 0.55M o (accounting for unresolved 
binaries gives F = -0.35). Thus, there are indica tions that 
the Class II IMFs are better describ ed by the IScalol Jl998h or 
the iRroupa. Tout & Gilmorel dl993h than by the more recent 
iKroupalfeOOll) . 



6. Cloud structure deduced from 1.3 mm and 
3.6 cm observations 

Figure ITU shows the IRAM 1.3 mm continuum mapping of 
a major part of the ISOCAM field. Four strong continuum 
sources are labelled in the figure. One is a multiple of three 
and another is elongated - all listed with coordinates and fluxes 
in Table The protosta r clustering scale of 0.12 pc found in 
the Serpens Cloud Core dKaas et all200 4^ is drawn in the map 
for size reference. Also shown are the positions of YSO can- 
didates found in Sect.0] Class I sources, flat-spectrum sources, 
and Class II sources. 

The contour plot of the 1.3 mm continuum emissio n out- 
lines t he position of the two ammonia cores found by Clark 
dl99ll) on each side of the optically visible stellar group 
Ser/G3-G6 quite well. MMS4 is a starless core coinciding with 
an absorption core seen in both ISOCAM filters. It is a candi- 
date prestellar core similar to the cores s tudied in absorption 
with ISOCAM bv iBacmann et alJ feOOOl) . MMS1 is located 
in the Ser/G3-G6SW ammonia core and is a multiple source 



Fig. 14. NOTCam 2.087 /im continuum image of Ser/G3- 
G6NE. The same field as Fig.lHl 



whose brightest components we have denoted as a, b, and c. 
Both MMS2, which is extended, and MMS3 are located in the 
Ser/G3-G6NE ammonia core. 

Morphologically, the 1.3 mm map is quit e similar to th e 
IRAM map of the Serpens Cloud Core dKaas et all l2004 . 
showing two separate clumps of which one has an elongated 
tail in the direction of the alignment of the clumps along an 
elongated density structure or a cloud filament. Here the elon- 
gated density structure oriented NE-SW is also well traced in 
the /f-band image of Fig.|2]as a pronounced drop in the surface 
density of faint stars. In the case of the Serpen s Cloud Core , 
the NW-SE oriented ridge was described in, e.g.. lKaasl(ll9*9 9a). 
Each clump contains a sub-structure. This region differs from 
the Cloud Core, however, in that the number of protostars is 
lower, and the 1.3 mm emission is in general less intense. The 
immediate interpretation is that we see no such intense burst 
of star formation here as in the Cloud Core, which is also in 
agreement with the Class I/Class II number fraction discussed 
at the end of Sect. fO 

6.1. Individual 1.3 mm and 3.6 cm sources 

The fluxes and positions of the IRAM 1 .3 mm and VLA 3.6 cm 
sources are given in Tables and [S] Taking into account the 
uncertainties in the ISOCAM positions (estimated in Sect. 12. 2> 
and the IRAM positions (~ 3" on the average and 5" at most), 
here we compare the mm and cm source positions with those 
of ISOCAM. 
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Table 7. The bright IRAM 1.3 mm sources with J2000 posi- 
tions good to about +3 arcsec. 



1^ 3.1X16 


^2000 


A™™ 

02000 


Peak S, , 
redK O l.3 m m 


nil o l.3mm 








(mjy/beam) 


(mJy) 


MMSl-a 


18 28 54.1 


00 29 29 


234. 


830 


MMSl-b 


18 28 54.3 


00 29 02 


153. 


700 


MMSl-c 


18 28 55.8 


00 29 45 


161. 


780 


MMS2 2 


18 29 06.5 


00 30 41 


404. 


1650 


MMS3 


18 29 09.2 


00 31 32 


399. 


1020 


MMS4 


18 29 31.2 


00 26 50 


94. 


420 



1 Integrated over an aperture with 37" diameter. 

2 This source is elongated. 



ISO-62 is probably the IR counterpart of MMSl-a, as the 
coordinates agree to 1.8". There is no A"-band detection, but 
L' = 1 1.4 mag, which gives a mid-IR SED index a A 1R 14 = In- 
consistent with a Class I source classification. The radio source 
VLA-#2 is positioned only 3.8" away from the mm source. The 
MMSl-b position differs by 8.5" from the IR coordinates of 
ISO-64, which is also a Class I source. Since this is just below 
the 3cr in positional uncertainty, it can not be excluded that 
ISO-64 is the infrared counterpart of MMSl-b. At the position 
of MMS 1 -c we find no IR counterpart, although we note that 
this region suffered strongly from persistency effects in both 
ISOCAM maps. There is, however, a 3.3 cr radio source (VLA- 
#3) only 1.4" to the south of MMSl-c. We tentatively suggest 
MMSl-c as a possible Class candidate. 

MMS2 is elongated, and its position determined by the 
peak emission is found between the ISOCAM positions of the 
two Class I sources ISO-90 and ISO-94. The MMS2 elonga- 
tion is roughly in the direction towards the mid-IR sources, and 
we interpret these to be the IR counterparts of the elongated 
MMS2 source, possibly an unresolved double. (The offsets be- 
tween the peak MMS2 coordinate and the ISOCAM positions 
are 7.8" and 12" from ISO-90 and ISO-94, respectively.) We 
also note that a radio source (VLA-#6) is found at the position 
of ISO-94 offset by only 6.2". We note that the Arnica/NOT po- 
sition of the Ks = 15.5 mag source (NOT-343) agrees within 
0.9" with the VLA position. 

MMS3 is separated from the faint mid-IR source ISO- 
101 by 6". ISO- 101 has no K-band detection, and in the L' - 
band only a very faint extended signal was seen. Judging from 
Fig.^]the position of ISO-101 seems to coincide with the ex- 
tended H 2 line emission. Also, the ISOCAM colour [14.3/6.7] 
is rather blue compared to the lower limit on the a 2 IR lA in- 
dex. This is likely explained by H2 line emission in the 6.7 
pm band and suggests that ISO-101 is a just a hot spot in the 
bright jet-like feature seen in the deep Ho line image in Fig. 1131 
rather than being a continuum source. Most probably, ISO-101 
is not the mid-IR counterpart of MMS3. There is a radio source 
(VLA-#7) less than 3" from the IRAM position of MMS 3. 
Based on this we propose that MMS3 is a new Class can- 
didate and probably the driving source of the bright outflow 
seen in the deep Hi line image in Fig. ^](see discussion in 
Sect.lO. 



Table 8. The VLA 3.6 cm sources with J2000 positions good to 
±0.5 arcsec. All except one are detections above 5<x. The beam 
size is given in Sect. 12.71 



V LA 


tt2000 


r 

02000 


c 

^3. 6cm 


Associations 


l#) 






(mJy/beam) 




1 


18 28 53.49 


00 29 11.9 


0.08 




2 


18 28 54.05 


00 29 32.7 


0.08 


MMSl-a, ISO-62 


(3) 1 


18 28 55.80 


00 29 43.6 


0.05 


MMSl-c 


4 


18 29 01.80 


00 29 54.4 


3.30 


ISO-82, Ser-G3 


5 


18 29 03.31 


00 30 57.9 


0.09 




6 


18 29 06.70 


00 30 33.6 


0.32 


ISO-94, MMS2 


7 


18 29 09.03 


00 31 30.6 


0.08 


MMS 3 


8 


18 29 11.14 


00 29 24.8 


0.12 


ISO- 107 


9 


18 29 25.52 


00 25 04.0 


0.14 





1 This is a 3.3<x detection only, and therefore only tentative. 



MMS4 has no IR sources in its vicinity, and as described 
above it is a candidate prestellar core. The VLA source #4 is 
a very bright radio source coinciding in position with ISO-82 
(CoKu Ser-G3), a Class II source that i s seen optically and wa s 
classified as a Classical TTauri star by Cohen & Kuhil dl979l) . 
The position of another radio source (VLA-#8) is only 4.2" 
away from ISO-107, a bright (K = 7.8 mag) star in the cen- 
tral part of the cluster that has no infrared excess, but could 
possibly be a YSO of type Class III. Three VLA sources (# 1, 
5, and 9) are without mm or IR counterparts and are probably 
extragalactic sources. 

Radio continuum emission at 3.6 cm has been detected for 
a number of YSO s - most of them in the Class or Class I 
protostellar phase (Eiro a et all2005h . Spectral indices obtained 
at cm wavelengths show that thermal free-free emission from 
ionised gas usually dominates, and it has been suggested that 
the radio emission from protostars is due t o thermal radio jets 
ionized by shocks in strong stellar winds JAnglada et all 1 998: 
Rei ourth et al. 2004). Although we have no information on the 
spectral index in the radio, this is probably the most likely in- 
terpretation for our VLA sources #2, 3, 6, and 7, which all 
correspond to Class I or Class sources, and of which two are 
associated with H2 outflows (se e next section). Fo r compari- 
son, in the Serpens Cloud Core lEiroa et al ] J2005I) found that 
as many as 10 of 16 YSOs with radio emission were in the 
Class to Class I phase, while only about 25% were as evolved 
as Class II sources. For this last group it has been suggested that 
non-thermal gyrosynchrotron emission from coronally active 
stars may explain the radio emission llSmith. Gtidel. & Benzl 
1999). We speculate that this could be the case for VLA # 4 and 
8, although we have no X-ray information for these sources. 
Determination of the radio spectral indices is needed to con- 
firm the above suggestions. 

6.2. Outflows traced by the 2. 122 11m S(1) line of H 2 

The Ser/G3-G6NE ammonia core was imaged with NOTCam 
in the 2.122 /im S(l) line of H2 (see Fig. II 3I >. Numerous 
Herbig-Haro (HH) objects with pure H2 line emission are 
found by comparison with the 2.087 //m continuum image in 
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iso-59 

Fig. 15. ALFOSC/NOT image of 900s integration through the 
H a filter (18 nm wide) in June 2002. Apart from some extended 
emission around the Ser/G3-G6 group, the only pure line emis- 
sion object seen in the optical is the already known HH 476. 
FOV=6.5'. 




Fig. 16. NOTCam H 2 line image (2.122 //m) of the Ser/G3- 
G6NE core. ISOCAM 6.7 fim contours are overlaid, and 
ISOCAM sources are labelled with the ISO ID number of 
Tabled 

Fig-El The morphology is complex, and proper motion of the 
knots is not yet available, but the features probably arise from at 
least two different bipolar collimated flows crossing each other 
along the line of sight. From the optical H a image in Fig. 1151 
it is evident that the flows are deeply embedded. We tentatively 
suggest that the driving sources are the two dense mm cores 
MMS2 and MMS3, the last one proposed to be a new Class 
candidate. 



Among the two bipolar flows traced by the numerous HH 
objects seen in the near-IR in Fig.^J one is a very long colli- 
mated flow oriented NE-SW and extending over the full FOV 
of the image. It probably arises from MMS2, either of the pos- 
sible multiple or as a consequence of disintegrating multiples, 
a mec hanism suggested by Reirjurth (2000) and Reip urth et alJ 
(2004). Its extension is at least 90" to both sides of the driv- 
ing source. We notice that the flow is slightly S-shaped. In 
the NE it ends in a bow-shock-like morphology (upper left 
of image). The optical Herbig-Haro object HH 476 that falls 
just outside the H2 image (but is shown in the larger field Hq, 
image in Fig. \15) may be related to this flow. HH 476 was 
detected by Ziener & Eisloffel ( 1999) in an optical [SII] sur- 
vey, and IWu etalJ ikOol" suggested its energy source to be 
IRAS 18265 +0028, which ISOCAM resolved in the two very 
red sources ISO-90 and ISO-94. Follow-up imaging is needed 
to determine the proper motion of the knots and to confirm our 
tentative interpretation. 

The remaining H2 emission features, i.e., the bright knots 
to the south of MMS3 and grossly aligned almost N-S, and also 
the faint extended emission seen to the north of MMS3 around 
the position of ISO- 106, most likely arise from another out- 
flow, one that is probably driven by MMS3, a new Class can- 
didate. The morphology is more complex with ringlike shapes, 
and more data is needed to interpret these structures. 

The only extended continuum radiation seen in Fig. 1141 is 
some faint nebulous emission in the vicinity of the two sug- 
gested driving sources, MMS2 and MMS3, and this is probably 
scattered light. 

7. Summary and conclusions 

An embedded cluster of YSOs was found from an ISOCAM 
survey around the optical stellar group Ser/G3-G6. Mid-IR 
photometry in two broadband filters centred on 6.7 and 14.3 
jum obtained with ISOCAM was combined with data from 
various telescope/instrument configurations to present a multi- 
wavelength study. This reveals low-mass YSOs in the evolu- 
tionary stages from Class to Class III. We propose 2 Class 
candidates, 5 Class I candidates, 5 flat-spectrum sources, 31 
Class II objects, and 3 Class Ills (with some IR-excess). Our 
selection criterium is IR excess, which means that our sam- 
ple is severly incomplete for Class III sources in general. This 
also means that we cannot estimate the disk fraction in the YSO 
population. The number fraction of Class I to Class II sources in 
the Ser/G3-G6 IR cluster (5/31) is quite typical for star-forming 
regions and comparable to the p Ophiuchi and Chamaeleon I 
clusters, i.e., quite different from the unusu ally large fractio n 
of Class Is found in the Serpens Cloud Core jKaas et all2004l) . 

The IRAM 1.3 mm continuum mapping detects several 
bright continuum sources in the two ammonia cores Ser/G3- 
G6NE and Ser/G3-G6SW, and one outside that is suggested as 
a candidate prestellar core. MMSl-c and MMS3 are proposed 
to be Class candidates. Both of them are also VLA 3.6 cm 
radio continuum sources. Deep NOTCam imaging in the 2.122 
jum line of H2 reveals signs of two complex outflows in the 
Ser/G3-G6NE core. A deep H a image shows that the flows are 
embedded and only the previously known HH 476 is seen opti- 
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cally. We suggest that the driving sources of these two outflows 
are MMS2 and MMS3. 

Comparing with existing pre-main sequence evolutionary 
models, we find that for any reasonable assumptions on age for 
the Class II sources, our sample extends well into the BD mass 
regime. The Class II sample was joined with the Class II sample 
from the Serpens Cloud Core to improve statistics. Model LFs 
calculated for a set of ages and I MFs give a best age of 2 Myr 
for which both the lScalol i 19981) and lKroupal fcOOll) IMFs are 
compatible with the observed LF down to the completeness 
limit at 0.08 L . The Class II LF can be described by co-eval 
star formation ~ 2 Myr ago, continuous star formation over the 
last 2 Myr, or a burst of star formation of various duration that 
took place 2 Myr ago. There is a peak in the observed LF that 
cannot be used to constrain the age because it coincides with 
the completeness limit. Nevertheless, the general shape of the 
LF does not permit a much younger or older population. The 
completeness limit must be shifted to 0.04 L Q to possibly dis- 
tinguish between different IMFs and star formation scenarios. 
This should be feasible with the deeper mid-IR surveys now 
available from Spitzer. 

Acknowledgements. We thank the referee, Dr Paul Harvey, for sug- 
gestions that led to substantial improvements of the paper. We ac- 
knowledge Dr. F.O. Clark for presenting his unpublished NH 3 map- 
ping to us before planning the ISOCAM observations. The ISOCAM 
data presented in this paper was reduced using "CIA", a joint develop- 
ment by the ESA Astrophysics Division and the ISOCAM Consortium 
led by the ISOCAM PI, C. Cesarsky, Direction des Sciences de 
la Matiere, C.E.A., France. This work is partly based on obser- 
vations made with the NOT, operated on the island of La Palma 
jointly by Denmark, Finland, Iceland, Norway, and Sweden, in the 
Spanish Observatorio del Roque de los Muchachos of the Instituto 
de Astrofisica de Canarias. One image presented here has been 
taken using ALFOSC, which is owned by the Instituto de Astrofisica 
de Andalucia (IAA) and operated at the NOT under an agreement 
between IAA and the NBIfAFG of the Astronomical Observatory 
of Copenhagen. This publication makes use of data products from 
the Two Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. This research has made use of the SIMBAD database, 
operated at CDS, Strasbourg, France, as well as SAOImage DS9, de- 
veloped by the Smithsonian Astrophysical Observatory. 



References 

Adams, F.C., Lada, C.J., & Shu, FH. 1987, ApJ 312, 788 
Anglada, G., Villuendas, E., Estalella, R., et al. 1998, AJ 116, 
2953 

Andre, P., Ward-Thompson, D., & Barsony, M. 1993 ApJ 406, 
122 

Andre, P., Ward-Thompson, D., & Barsony, M. 2000, in 
Protostars and Planets IV, ed. V. Mannings, A. P. Boss, & S.S. 
Russell (Tucson: University of Arizona Press), p. 59 
Andre, P. & Montmerle, T. 1994, ApJ 420, 837 
Bacmann, A., Andre, P., Puget, J.-L., et al. 2000, A&A 361, 
555 



Baraffe, I., Chabrier, G., Allard, F, & Hauschildt, PH. 1998, 

A&A 337, 403 
Bessell, M.S., & Brett, J.M. 1988, PASP 100, 1134 
Blommaert, J., Metcalfe, L., Altieri, B., et al. 2000, 

Experimental Astronomy, 10, 241 
Bonnell, I.A., Larson, R.B., & Zinnecker, H. 2006, in 

Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. 

Keil, in press, astro-ph/0603447 
Bontemps, S., Andre, P., Kaas, A.A., et al. 2001, A&A 372, 

173 

Broguiere, D., Neri, R., & Sievers, A. 1995, NIC bolometer 

users guide (IRAM internal report) 
Casali, M.M., & Hawarden, T.G 1992, UKIRT Newsletter 4, 

33 

Clark, F.O. 1990, private communication 
Clark, F.O. 1991, ApJS 75, 611 

Cesarsky C.J., Abergel, A., Agnese, P., et al. 1996, A&A 315, 
L32 

Cohen, M., & Kuhi, L.V. 1979, ApJS 41, 743 
Curiel, S., Rodriguez, L.F, Moran, J.M., & Canto, J. 1993, ApJ 
415, 191 

Cutri, R.M., Skrutskie, M.F, Van Dyk, S., et al. 2003, 
Explanatory Supplement to the 2MASS All Sky Data 
Release, 

http://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html 
Dame, T.M., & Thaddeus, P. 1985, ApJ 297, 751 
D'Antona F, & Mazzitelli I. 1998, in Brown Dwarfs and 

Extrasolar Planets, ed. R. Rebolo, E.L. Martin, & M.R. 

Zapatero Osorio, ASP Conf. Proc. 134, p. 442 
Emerson, D.T., Klein, U., & Haslam, C.G.T. 1979, A&A, 76, 

92 

Eiroa, C, Torrelles, J.M., Curiel, S., & Djupvik, A.A. 2005, AJ 
130,643 

Eiroa, C, Djupvik, A.A., & Casali, M.M. 2006, in "Handbook 

of star formation regions", ed. B. Reipurth, in press 
Folha, D.F.M., & Emerson, J.P 1999, A&A 352, 517 
Froebrich, D., Ray, T.P, Murphy, G.C., & Scholz, A. 2004 
A&A432,L67 

Greene, T.P, Wilking, B.A., Andre, P., Young, E.T., & Lada, 

C.J. 1994, ApJ 434, 614 
Griffin, M.J., & Orton, GS. 1993, Icarus, 105, 337 
Harvey, P.M., Chapman, N., Lai, S.-P, et al. 2006, ApJ, in press 
Hunt, L., Mannucci, F, Testi, L., et al. 1998, AJ 1 15, 2594 
Kaas, A.A. 1999a, AJ, 118, 558 

Kaas, A.A. 1999b, PhD Thesis, Stockholm University, 
Akademitryck AB, Edsbruk, ISBN 91-7153-923-9 

Kaas, A.A., Olofsson, G, Bontemps, S., et al. 1999, in The 
Universe as seen by ISO, ed. P. Cox & M.F. Kessler, ESA 
SP-427, p. 493 

Kaas, A.A., & Bontemps, S. 2000, in From Darkness to 
Light: Origin and Evolution of Young Stellar Clusters, ed. 
T. Montmerle & P. Andre, ASP Conf. Proc. 243, p. 367 

Kaas, A.A., Olofsson, G, Bontemps, S., et al. 2004, A&A 421, 
623 

Kessler, M.F, Steinz, J.A., Anderegg, et al. 1996, A&A, 315, 
L27 

Koornneef, J. 1983, A&A, 128, 84 



Djupvik et al.: The young embedded cluster around Serpens/G3-G6 



Kreysa, E., Gemiind, H.-R, Gromke, et al. 1998, in Advanced 
Technology MMW, Radio, and Terahertz Telescopes, Ed. 
T.G. Phillips, Rroc. SPIE, 3357, 319 

Kroupa, P., Tout, C, & Gilmore, GF. 1993, MNRAS 262, 545 

Kroupa, P. 2001, MNRAS 322, 231 

Lada, C.J. 1987, in IAU Symp. 1 15, Star Forming Regions, ed. 

M. Peimbert & J. Jugaku (Dordrecht: D. Reidel Publishing 

Co.), p. 1-17 
Lada, C.J., & Wilking, B.A. 1984, ApJ 287, 610 
Lawson, W.A., Crause, L.A., Mamajek, E.E., & Feigelson, 

E.D. 2002, MNRAS 329, L29 
Lyo, A.-Ran, Lawson, W.A., Mamajek, et al. 2003, MNRAS 

338,616 

Meyer, M.R., Calvet, N., & Hillenbrand, LA. 1997, ApJ 114, 
288 

Olofsson, G., Kaas, A.A., Nordh, L., et al. 1998, in Brown 
Dwarfs and Extrasolar Planets, ed. R. Rebolo, E.L. Martin, 
& M.R. Zapatero Osorio, ASP Conf. Proc. 134, p. 81 

Olofsson, G, Floren, H.-G., & Lagage, P.O. 1999, Imaging in 
the Thermal Infrared, in Astrophysics with the NOT, ed. H. 
Karttunen, & V. Piirola (University of Turku), p. 37 

Persi, P., Palagi, F, & Felli, M. 1994, A&A 291, 577 

Persi, P., Marenzi, A.R., Olofsson, G., et al. 2000, A&A 357, 
219 

Reipurth, B. 2000, AJ 120, 3177 

Reipurth, B., Rodriguez, L.F, Anglada, & G., Bally, J. 2004, 

AJ 127, 1736 
Salpeter,E. 1955, ApJ 121, 161 

Scalo, J. 1998, in The Stellar Initial Mass Function, ed. G. 
Gilmore, D. Howell, Astronomical Society of the Pacific vol 
142, p. 201 

Smith, K., Gudel, & M., Benz, A.O. 1999, A&A 349, 475 
Straizys, V., Cernis, K, & Bartasiute, S. 2003, A&A 405, 585 
Strom, K.M., Strom, S.E., & Merrill, K.M. 1993, ApJ 412, 233 
van der Bliek, N, Olofsson, G, Floren, H.-G, et al. 2000, 
SIRCA, The Stockholm Infrared Camera, in The NOT in the 
2000's, ed. N. Bergvall, L. Takalo, & V. Piirola (University 
of Turku), p. 151 
Whittet, D.C.B. 1988, The observed properties of interstellar 
dust in the infrared, in Dust in the Universe, ed. M.E. Bailey 
& DA. Williams (New York: Cambridge University Press), 
p. 25-53 

Wu, J.W., Wu, Y.F, Wang, J.Z., & Cai, K. 2002, Chin. J. 

Astron. Astrophys. Vol 2, No 1, 33-42 
Zhang C.Y., Laureijs R.J., Clark F.O., & Wesselius PR. 1988, 

A&A 199, 170 
Ziener, R., & Eisloffel, J. 1999, A&A 347, 565 
Zinnecker, H., McCaughrean, M.J., & Wilking, B.A. 1993, 

in Protostars and Planets III, ed. E. Levy, & J.I. Lunine, 

(Tucson: University of Arizona Press), p. 429 



